Dynamical pinning and non-Hermitian mode transmutation in the Burgers equation 
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We discuss the mode spectrum in both the deterministic and noisy Burgers equations in one 
dimension. Similar to recent investigations of vortex depinning in superconductors, the spectrum is 
given by a non-Hermitian eigenvalue problem which is related to a 'quantum' problem by a complex 
gauge transformation. The soliton profile in the Burgers equation serves as a complex gauge field 
engendering a mode transmutation of diffusive modes into propagating modes and giving rise to a 
dynamical pinning of localized modes about the solitons. 
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The noisy Burgers equation and the related Kardar- 
Parisi-Zhang (KPZ) equation provide a continuum de- 
scription of an intrinsically nonequilibrium noise-driven 
system. As such they delimit an interesting class of sys- 
tems far from equilibrium. Specifically the equations ap- 
ply to the growth of an interface either due to a random 
drive or subject to random environments. 

In the case of one spatial dimension, which is our 
concern here, the Burgers equation for the local slope 
U = V/i has the form |Lj 
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vV 2 u + Vr] , 



(1) 



The related KPZ equation for the height h is dh/dt = 
i/V 2 h + (A/2)(Vft.) 2 +rj. The damping constant v charac- 
terizes the linear diffusive term. The coupling strength A 
controls the nonlinear growth or mode coupling term. 
The noise r\ is spatially short-ranged Gaussian white 
noise correlated according to (ri(xt)?] (00)) = A5(x)S(t) 
and characterized by the strength A. 

The stochastic equation ([!]) and its KPZ version have 
been studied intensively in particular in recent years and 
much insight concerning the pattern formation and scal- 
ing properties engendered by these equations has been 
gained on the basis of i) field theoretical approaches || , 
ii) mapping to directed polymers JM , and iii) mapping to 
the asymmetric exclusion model B. 

In recent works [pi we advanced a Martin-Siggia-Rose 
based canonical phase space approach to the noisy Burg- 
ers equation (Q). This method applies in the weak noise 
limit A — ► and replace the stochastic Burgers equation 
(El) by two coupled deterministic mean field equations 



d_ 
at 
d_ 

dt 



AwV u= vV 2 u - V 2 p 



- AwV p= -vV z p , 



(2) 
(3) 



for the slope field u and a canonically conjugate noise 
field p, characterizing the noise r\. The field equations 
derive from a principle of least action with Hamiltonian 
density TL = p(vV 2 u + XuVu — (1/2) V 2 p) and determine 
orbits in a canonical phase space spanned by u and p. 
Moreover, the action associated with a finite time orbit 
from u' to u, S = J '", dtdx(pdu/dt — TC), provides di- 
rect access to the transition probability P{v! — ► u, t) oc 
exp[— S/A] and associated correlations; an important as- 
pect which was pursued in || . 

On the 'zero noise' manifold for p = the field equa- 
tion (H) reduces to the noiseless Burgers equation 



vV 2 u , 



(4) 



which as a nonlinear evolution equation exhibiting tran- 
sient pattern formation has been used to model 'turbu- 
lence' and for example galaxy formation || , see also || . 

In this Letter we discuss two new features associated 
with the pattern formation in the Burgers equation in 
both the noiseless case (^) and the noisy case in terms of 
(||) and (H) . We focus on the interplay between localized 
nonlinear soliton modes and superposed linear modes and 
show the soliton-induced i) mode transmutation of dif- 
fusive modes into propagating modes and ii) dynamical 
pinning of linear modes about the solitons. Details will 
appear elsewhere. 

It is a feature of the nonlinear growth terms that the 
field equations (§) and (||) admit nonlinear localized soli- 
ton solutions, in the static case of the kink-like form 



= jm tanh[fc s 



k s — 



Xu 
2^' 



±1 . 



(5) 



The index \i labels the right hand soliton for (i = 1 with 
p s = 0, also a solution of the damped noiseless Burgers 
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equation for rj = 0; and the noise-excited left hand soli- 
ton for /J, = — 1 with p s = 2vu a , a solution of the grow- 
ing (unstable) noiseless Burgers equation for v — > — v. 
The amplitude-dependent wavenumber fc s sets the in- 
verse soliton length scale. Noting that the field equa- 
tions (||) and (|§|) are invariant under the slope-dependent 
Galilean transformation 



Xuot , 



(6) 



propagating solitons are generated by the Galilean boost 
(^J). Denoting the right and left boundary values by w+ 
and it_, respectively, the propagating velocity is given by 
the soliton condition 



u + + u- = —2v/\ 



(7) 



It follows from the quasi-particle representation ad- 
vanced in ||, see also p0| , that a general interface slope 
profile u — u s + Su at a particular instant can be rep- 
resented by a dilute gas of solitons amplitude-matched 
according to (^|) with superposed linear modes Su. For a 
configuration consisting of n solitons we then have 



2v n 

u s = — k p tanh — v p t — x p ) , 

P =i 



(8) 



where we have introduced the mean amplitude of the 
p-th soliton k p = (A/4z/)(u p +i — u p ), u p+ i and u v are 
the boundary values. The velocity of the p-th soli- 
ton is v p = — (A/2)(u p+ i + Up), x p is the center of 
mass, and we are assuming vanishing boundary condi- 
tions i*i = u n+ \ — 0. Note that the configuration (|^) 
is only valid at times in between soliton collisions; the 
interface changes dynamically subject to the conserva- 
tion of energy, momentum, and total area J dx u under 
the slope profile. The number of solitons, however, is 
not conserved, see HJh}. In Fig 1. we have depicted an 
n-soliton configuration 

Us 



p+1 



Jn+1 



U n 



FIG. 1. We depict an n-soliton slope configuration of 
a growing interface. The p-th soliton moves with velocity 
v p — — (A/2)(u p +i +u p ), has boundary value u+ and u_, and 
is centered at x p . The arrows on the horizontal inter-soliton 
segments indicate the propagation of linear modes. 

In order to discuss the linear mode spectrum about the 
soliton configuration u s it is convenient to introduce the 
shifted noise field ip 



p = v(u — if) 



(9) 



The field equations (§) and (0) then assume the symmet- 
rical form, also discussed in , 
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m 



AuVj u = v\J 2 p , 
AuV I ip — v\7 2 u . 



(10) 
(11) 



In the linear Edwards- Wilkinson case [[11J for A = the 
field equations readily support a diffusive mode spec- 
trum of extended growing and decaying modes, u ± 
p oc exp(=pi^fc 2 t) exp(ikx), i.e., u = [Aexp(— vk 2 t) + 
B exp(uk 2 t)] exp(z/cx), consistent with the phase space 
interpretation discussed in Q . 

Expanding about u s and the associated noise field (p s 
= M< for fi = ±1), 



2v 

— \k p \ tanh \ k p 

P =i 



Vr,t ■ 



(12) 



u = u s + 8u and tp — p s + Sip, the superposed lin- 
ear mode spectrum is governed by the coupled non- 
Hcrmitian eigenvalue equations 



- Au s V ) Su = v\J 2 8p + X(Vu s )6u , (13) 



d_ 

dt 

\ 

— - Xu s \7 j Sip = vX7 2 5u + X(Vp s )Su . (14) 



In the inter-soliton matching regions of constant slope 
field Vu s = V(f s — and the equations (|l^) and 
( |l4] ) decouple as in the Edwards- Wilkinson case. Set- 
ting u s — u and and searching for solutions of the 
form Su, Sip oc exp(-Ekt) exp(ikx) we obtain Su ± 
Sip oc exp(— E^t) exp(ikx), i.e., Su = [Aexp(—E£t) + 
B exp(— EZt)] exp(ikx), where the complex spectrum 
characteristic of a non-Hermitian eigenvalue problem is 
given by 



E. 



±isk — iXuk 



(15) 



Introducing the phase velocity v = Xu the Su mode cor- 
responds to a propagating wave with both a growing and 
decaying component 



Su oc (Ae 



-vh t 



uh t\ ik(x+vt) 



Be vk *)e 



(16) 



The presence of the nonlinear soliton profile thus gives 
rise to a mode transmutation in the sense that the diffu- 
sive mode in the Edwards- Wilkinson case is transmuted 
to a propagating mode ( |l6| ) in the Burgers case. As in- 
dicated in Fig. 1 the linear mode propagates to the left 
for u > and to the right for u < 0. We note in par- 
ticular that for a static right hand soliton u± = ±it and 
the mode propagates towards the soliton center which 
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thus acts like a 'sink'; for a static left hand soliton the 
situation is reversed, the mode propagates away from the 
soliton which in this case plays the role of a 'source'. 

In the soliton region the slope field varies over a scale 
give n by k^ 1 and we must address the equations ( [l3"|) and 
(|l4|). Introducing the auxiliary variables 



SX ± = Su±Sip , 

they take the form 

dSX ± 
dt~ ~ 

±D (±^ s ) SX* - ^(Vm s ± V(p s )6XT , 



(17) 



(18) 



where D(±\u s /2v) is the 'gauged' Schrodinger operator 



2v 



2v 4i/ 



2 f?<P. (19) 



for the motion of a particle in the potential {\/Av)u 2 — 
(\/2)\7ip s subject to a gauge field {\/2v)u s given by u s . 

In the regions of constant slope field Vit s = V</? s = 
0, u s = u, D{±\u s /2v) -> -i/(V ± \uj2v) 2 + 
(\/Av)u 2 and searching for solutions of the form SX^ cx 
cxp(— Ekt) exp(«/cx) we recover the spectrum (|l5|). In 
the soliton region V</?^ = /iVm^, fi = ±1, and one of 
the equations (^) decouples driving the other equation 
parametrically. 

In order to pursue the analysis of ( |l8| ) we note that 
the gauge field \u s j2v in the Schrodinger operator can 
be absorbed by means of the gauge or Cole-Hopf |l2] ] 
transformation 



U = exp 









f dx u s 


2v J 





(20) 



Using the relation D(±\u s /2v) = U^D^U^ 1 we ob- 
tain the Hermitian eigenvalue equations 



dSX ± 

±U ±1 D(0)U Tl 6X ± - -(Vu a ± Vtps)5XT 



(21) 



which are readily analyzed in terms of the spectrum of 
D(0) discussed in @. 

The exponent or generator in the gauge transforma- 
tion ( |20[ ) samples the area under the slope profile u s up 
to the point x. For x — * oo U — > exp[— XM/2iA, where 
M = J dx u s is the total area; according to (jl|) or 
M is a conserved quantity. In terms of the height field 
h, u = Vh, M = h(+L) — h(—L) is equal to the height 
offset across a system of size L, i.e., a conserved quan- 
tity under growth. Inserting the soliton profile u s (Jsj) the 
transformation U factorizes in contributions from local 
solitons, i.e., 



U = II U* igIik » , U p = cosh" 1 fc p (a 



v n t 



(22) 



Focusing on a particular soliton contribution to the inter- 
face with boundary values it+ and u_ and for convenience 
located at x p = the analysis is most easily organized 
by first performing a Galilean transformation (|6[) to a lo- 
cal rest frame by shifting the slope field by (u+ + w_)/2, 
corresponding to the velocity given by (0). The static 
soliton is then given by (^|), u± — ±u, and 



D(0) = -vV 2 + vk 2 s [l - 2/ cosh 2 k s x] 



(23) 



describing the motion of a particle in the attractive po- 
tential — 2^fc 2 /cosh 2 k s x whose spectrum is known. 

Denoting the eigenvalue problem D(Q)^ n — iln^ n the 
spectrum of D(0) is composed of a zero-energy CIq = lo- 
calized state cosh" 1 k s x, yielding the soliton trans- 
lation mode lifting the broken translational symmetry, 
and a band ^k °c exp(ik s x)sk(x) of extended phase- 
shifted scattering modes with energy 



n k = v{k 2 + k 2 ) . 



(24) 



Sk{x) = (k+ikg tanhk s x)/(k — ik s ) is a modulation of the 
plane wave state; for an oo st(i) -» exp(i<5fc), where 5k 
is the phase shift of the wave. 

Inserting Vtp% = /xVmJ the fluctuations 

<5X ± = U^SX ± , U p = cosh" 1 k s x (25) 

then satisfy the Hermitian eigenvalue equations 

dSX* 



dt 



±D(Q)SX ± - v{n ± 1)6X* , 



(26) 



which decouple and are analyzed by expanding SX^ on 
the eigenstates \& n . For the right hand soliton for /j, = +1 
and focusing on the plane wave component, ignoring 
phase shift effects, we obtain in particular the fluctua- 
tions 



SX+ = (Ae- Ukt + Be 

SX~ = B — e Ukt e lkx coshk s x 
v 



e"" x cosh 1 k s x , 



(27) 
(28) 



For real k the modes SX are diffusive and the spectrum 
ilk = v(k 2 + k 2 ) exhibits a gap vk 2 proportional to the 
soliton amplitude squared. Moreover, the gauge transfor- 
mation Up gives rise to a spatial modulation of the plane 
wave form which allows us to extend the spectrum by an 
analytical continuation in the wavenumber k. In particu- 
lar by setting k — ► k^f ik s for SX^ and noting that 5X^ 
decouple for i > i; s _1 we have fife — > vk 2 =p 2ivkk s and 
exp(ikx) cosh^ 1 k s — > const and we achieve a matching 
to the extended propagating modes in the inter-soliton 
regions. A similar analysis applies to the left hand soli- 
ton for fx = —1 and the gauge transformation ( p2| ) allows 
for a complete analysis of the linear fluctuation spectrum 
about the multi-soliton configuration u s . 
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The last issue we wish to address is the fluctuation 
spectrum in the noiseless case, extending the analysis in 
Pj . In this case we only have the right hand soliton for 
// = +1 and the soliton and associated fluctuations lie on 
the zero- noise manifold for p — 0. There is no coupling 
to the 'noisy' modes i.e., u = ip, 8u — Sp — 5X + /2, 
and 8X~ = 0, and the fluctuations are given by ( gjj ) for 
B = (ignoring phase shift effects) 



Su = 5X+/2 oc e 



^e ikx cosh -1 fc, 



(29) 



It is an essential feature of the non-Hermitian eigen- 
value problem ( |l8| ) characterizing noisy nonequilibrium 
growth, and in the noiseless case of transient growth 
yielding (29) that the real spectrum (|24|) can be extended 
into the complex eigenvalue plane. This is due to the en- 
velope cosh -1 k s x which gives rise to a spatial fall off. 
Setting k — > k + in, flk — * Ef. ;K , where 



v{k 2 



2ivkn 



(30) 



For |«| < k s we have a band of localized fluctuations 
dynamically pinned to the soliton. The modes are ex- 
ponentially damped with a damping constant given by 
the real part of Ej. tK , KeE^^ = v(k 2 — k 2 ). The imag- 
inary part of Ek,Ki Irf-E'fc.K — 2vkn, combined with the 
phase ikx yields a propagating wave with phase veloc- 
ity 2vk, finally the spatial range of the mode is given by 
(k s — k) . For k = the spectrum is real, the phase ve- 
locity vanishes, the range is fcj 1 , and the localized mode 
is symmetric and purely diffusive with a gap k 2 . For 
k = k g , the borderline case, the fluctuations are extended 
in space and purely propagating with a gapless spectrum 
vk 2 . For intermediate n values the modes are propagat- 
ing with localized envelopes. In Fig. 2 we have depicted 
the complex eigenvalue spectrum. 

Im E k:K 
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k,K 



FIG. 2. We depict the the complex eigenvalue spectrum 
for the damped modes in the noiseless Burgers equation. On 
the boundaries k, — ±fc s the modes are extended and prop- 
agating. The shaded area indicates propagating modes with 
localized envelopes. For k = the mode is symmetrical and 
purely diffusive. 

In this Letter we have analyzed two aspects of the 
interplay between linear superposed modes and nonlin- 
ear soliton excitations in the noiseless and noisy Burgers 



equations describing transient and stationary nonequilib- 
rium growth, respectively. Both aspects are intimately 
related to the non-Hermitian character of the eigenvalue 
problem. The first aspect is a linear mode transmutation 
where the diffusive non-propagating modes in the linear 
Edwards- Wilkinson case owing to the solitons are trans- 
muted to propagating extended modes in the nonlinear 
Burgers case. The second aspect is a dynamical pinning 
of a band of localized modes to the solitons. We finally 
note that similar aspects are also encountered in recent 
work on the transverse Meissner effect and flux pinning 
in superconductors. Here the uniform gauge field is given 
by the transverse magnetic field wheras in our case the 
nonlinear soliton profile provides the nonuniform gauge 
field. 

Discussions with A. Svane, K. M0lmer, J. Hertz, B. 
Derrida, M. Lassig, J. Krug and G. Schiitz are gratefully 
acknowledged. 
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